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ABSTRACT: During the aerosolization process at the sea surface, halides are 
incorporated into aerosol droplets, where they may play an important role in 
tropospheric ozone chemistry. Although this process may significantly 
contribute to the formation of reactive gas phase molecular halogens, little is 
known about the environmental factors that control how halides selectively 
accumulate at the air— water interface. In this study, the production of sea spray 
aerosol is simulated using electrospray ionization (ESI) of 100 nM equimolar 
solutions of NaCl, NaBr, Nal, NaN0 2 , NaN0 3 , NaC10 4 , and NaI0 4 . The 
microdroplets generated are analyzed by mass spectrometry to study the 
comparative enrichment of anions (J x ~) and their correlation with ion 
properties. Although no correlation exists between f x ~ and the limiting 
equivalent ionic conductivity, the correlation coefficient of the linear fit with 
the size of the anions R x , dehydration free-energy A G dehyd , and polarizability 
a, follows the order: R x ~ 2 > R x - -1 > &x“ > AG dehyd > a. The same pure physical process is observed in H 2 0 and D 2 0. The factor 
f x ~ does not change with pH (6.8— 8.6), counterion (Li + , Na + , K + , and Cs + ) substitution effects, or solvent polarity changes in 
methanol— and ethanol— water mixtures (0 < jc H2 q < l). Sodium polysorbate 20 surfactant is used to modify the structure of the 
interface. Despite the observed enrichment of I - on the air— water interface of equimolar solutions, our results of seawater mimic 
samples agree with a model in which the interfacial composition is increasingly enriched in 1“ < Br - < Cl - over the oceanic 
boundary layer due to concentration effects in sea spray aerosol formation. 


Seawater Ratios y Equimolar Ratios 



■ INTRODUCTION 

The fractionation of species occurring at the air— water interface 
of the ocean is an intriguing problem and the focus of 
continuous experimental and theoretical studies that aim to 
understand the fundamental role of aerosols and their effect in 
radiative forcing. 1 7 The process of seawater aerosolization 
during bubble bursting generates negatively charged droplets 
(for diameters between 0.2 and 2.5 //m) 8 and establishes a 
negative current of ions into the atmospheric boundary 
layer. 9 ' 10 Many key processes related to the transfer of anions 
such as halides (e.g., sea salt formation) and activation of 
halogens such as Br 2 to participate in the photochemical 
depletion of tropospheric ozone remain poorly understood. 11 
The depletion of tropospheric 0 3 (g) molecules in the polar 
regions 2,13 has been linked to interfacial reactions that convert 
inert halide ions (e.g., Br - ) into molecular halogens (e.g., Br 2 
and I 2 ) and reactive halogen species, RHS (e.g., Br*, BrO*, and 
IO*). Once generated, RHS deplete ozone ([0 3 ] « 40—60 
ppbv) in the Arctic boundary layer to ~0 ppbv. 12,14-16 A 2-fold 
important unsolved question is how are salts transferred from 
the ocean and then oxidized to become RHS in the air. 11,17 

The air— water interface enables unique reaction pathways 
that might be altered due to the high surface area available as 
compared to the bulk. It has been proposed that the affinity of 
larger anions toward the air— water interface results from their 
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rejection by the medium via collective dispersive interac- 
tions. 18 ' 19 Recently, a new theory that takes into account all the 
moments of the ionic charge distribution, and not just the 
dipole, predicted that larger anions are adsorbed at the 
interface, while the metal cations are repelled from it. 20 
Simulations at higher concentrations can be interpreted in 
terms of ion-specific ion— surface and ion— ion excluded volume 
correlations.” 1 Computational methods coupled with exper- 
imental results showed that, for dilute bulk concentrations, soft, 
polarizable monovalent anions such as I - have a tendency to be 
present at the interface due to the dipole— induced dipole 
interaction. 22 ' 23 Recently, electrospray ionization (ESI) mass 
spectrometry (MS) in the negative mode was proven to be a 
useful technique to study the accommodation of anions at the 
interface, and chemical reactions occurring during the 
formation of aerosol mimics. 18,24-27 Studies for >100 pM 
solutions described the size dependent fractionation of the 
outermost layers of aerosolized droplets using an ESI probe, for 
ions with mass to charge detection limit m/z > 50 amu. 18,24 In 
consequence, smaller anions such as chloride and nitrite have 
not been studied so far. Several mechanisms have been 
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proposed to control the selective enrichment of anions at the 
air— water interface. The mechanisms are not independent and 
include factors such as (1) ion size, 20 ' 24 (2) ion polarizability,” 
and (3) ion dehydration free energy. 18,28 

In this article, we study the distribution of seven monovalent 
anions upon ESI aerosolization, for 100 nM equimolar 
solutions, what is 1000 times more diluted than previously 
tested, 18 '” 4 as well as in mixtures that mimic their actual 
composition in seawater.” 0- 33 The fractionation of the 
environmental relevant and small Cl - and N0 2 - ions is 
monitored for the first time simultaneously with Br - , N0 3 - , I - , 
C10 4 - , and I0 4 - . The process is studied at the interface of 
water, D 2 0, water— methanol, and water— ethanol mixtures, 
which are largely covered by —OH, — OD, — CH 3 , and 
— CH 2 CH 3 groups, respectively. 34-37 The use of nonionic 
polysorbate 20 surfactant provides direct physical insight to 
understand the effect of polarizability (a) above the critical 
micelle concentration (CMC). We show that anion fractiona- 
tion in micrometer size droplets depends on the actual 
concentration ratios of the species present in surface ocean 
waters. This work also suggests that the reciprocal square of ion 
radius can be used as a better predictor for the enrichment of 
anions and that the process at relevant salts molar ratio is still 
dependent on ion size. 

■ EXPERIMENTAL SECTION 

Chemicals. Sodium salts were chosen to study the 
fractionation of anions at the air— liquid interface during the 
aerosolization process: NaCl (99%), Nal (99.99%), and 
NaN0 2 (99%, all from Fisher Scientific); NaI0 4 (99%) and 
NaBr (99.5%, both from Acros Organics); and NaN0 3 (99%) 
and NaC10 4 (99.8%, both from MP Biomedicals). In addition, 
potassium, lithium, and cesium halide salts were used to study 
the effect of cation substitution during aerosolization: KC1 
(99%), Lil (99%), and CsCl (99.9%, all from MP Biomedical); 
LiCl (>98%) and LiBr (99%, both from Alfa Aesar), CsBr 
(99.9%) and Csl (99.9%, both from Acros Organics); KBr 
(99%, Aldrich); and KI (99%, BDH Chemicals). 

Preparation of Solutions. A 1.00 mM stock solution of 
each salt was freshly prepared daily by direct weighing and 
dissolution in a calibrated volumetric flask. The stock solutions 
were used to prepare 400 nM equimolar mixtures that 
underwent a 4 time dilution during infusion to a final 
concentration of 100 nM in each anion (unless indicated 
otherwise). Solutions were prepared in ultrapure water 
provided by an Elga Purelab flex (Veolia) water-polishing 
system (18.2 MO cm, <5 ppb organics (TOC), <1 CFU/mL), 
D 2 0 (99.9%, Cambridge Isotope), methanol (99.9%, Fisher 
Optima), and ethanol (99.5%, Sigma Aldrich). 

Equimolar solutions in water were used to study the 
dependence of fractionation on (l) concentration, (2) acidity 
(6.8 < pH < 8.6), (3) solvent polarity with alcohol mixtures (0 
< »h 2 o ^ l)j and (4) the presence of 10 -8 — 10 -3 M polysorbate 
20 (100%, Sigma Aldrich) as a nonionic surfactant. Samples at 
representative concentration of monovalent anions in surface 
seawater 29-33 were prepared. 

The pH of a 1 L solution in water, with 100 nmoles of each 
sodium halide studied, initially at 5.96, was adjusted by the 
addition of 1.000 X 10 -1 M NaOH (Acros Organics) and 
measured with a pH meter (Mettler Toledo) calibrated at pH 
4.01 and 7.00 or 7.00 and 10.01 with buffer solutions (Orion). 


ESI-MS Description. The process of ESI has been 
efficiently used for many years as a mechanism of aerosolization 
of nonvolatile molecules such as highly polar, ionic, and 
macromolecular analytes, into the gas-phase region of a mass 
spectrometer. 38 In our ESI system (Scheme l), the dissolved 


Scheme 1. Electrospray Ionization Mass Spectrometer (ESI- 
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ionic compounds start their journey to the surface in a liquid 
solution (H,0, D 2 0, H 2 0/CH 3 0H, or H 2 0/CH 3 CH 2 0H 
mixtures) that is nebulized into an atmospheric pressure ion 
source region. Pneumatically assisted electrospray converts the 
solution into small droplets by the combination of a strong 
electric potential (1.5 kV) between needle (a stainless steel 
capillary with 127 fim ID and 230 /(m OD) and counter 
electrode and a high-flow of N 2 nebulizing gas (P = 70 psi, flow 
= 12 L min -1 ). The capillary carrying the sample and solvent 
molecules is surrounded by two concentric tubes carrying the 
nebulizing gas; the first one (330 fim ID) assists in the 
aerosolization, and the second tube transports the sheath gas to 
help in the desolvation process of the liquid jet. The larger 
velocity of the N 2 nebulizing gas (^ gas = 3.12 X 10 4 cm s -1 ) over 
the liquid forming the plume of droplets (liquid = 26.3 cm s -1 ) 
contributes to the formation of fine aerosol droplets. The 
operation of the ESI probe in the negative-ion mode generates 
small negatively charged parent droplets with radius R d , which 
subsequently form a plume of solvated molecular ions that 
travel between the needle and the entrance cone that acts as a 
counter electrode. Neutral solvent molecules undergo evapo- 
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ration from the surface of the charged droplets, causing the 
droplets to shrink in size, shortening the distance of the 
negative charges at the surface of the newly formed smaller 
droplets. Thus, the ions are sputtered off the surface as the 
charge builds up, making the experiment surface sensitive. 24 
The first generation daughter droplets of smaller size with a 
radius R m experience charge crowding on the surface that 
generates repulsion forces that overcome the surface tension of 
the liquid. In consequence, the parent droplets disintegrate into 
daughter microdroplets that are 10 times smaller as indicated 
by the ratio 1^/Rj < 0.1. The ionized species are attracted 
toward the lower pressure orifice of the entrance cone set to 
—80 V, that as a counter electrode also attracts the ions. The 
ions travel into the first vacuum stage through the entrance 
capillary, where the majority of desolvation takes place and, 
finally, to a quadrupole mass analyzer (accurate to 0.17 amu), 
where they are detected at specific m/z ratios. The overall time 
from the formation of original droplets until ion detection is <1 
ms. The extracted ion count at the specific m/z of interest 
showed a direct correlation with concentration. However, 
converting the ion count into a fractionation ratio yields a 
concentration independent variable that reports the behavior of 
the ions on the interface rather than in the bulk of the 
solution. 18 '” 4 

The nebulizer voltage of the ESI probe was optimized to 1.5 
kV for maximum ion count, within a range that showed the 
same fractionation behavior for all the ions. The cone voltage, 
capable of producing collisional induced dissociation (CID) of 
the ions, was optimized to —80 V to obtain the highest ion 
count for all the ions of interest; although I0 4 - could be 
converted to I0 3 - at this voltage, its <5% loss is considered 
minor. In summary, the experimental conditions were (unless 
indicated otherwise in a figure) drying gas temperature, 250 °C; 
nebulizer voltage, —1.5 kV; cone voltage, —80 V; and nebulizer 
pressure, 70 psi. Mass spectra were acquired between 30 and 
200 amu, and the peaks at m/z 35 and 37 ( 3S,37 C1 _ ), 46 
(N0 2 “), 62 (N0 3 “), 79 and 81 ( 79,81 Br - ), 99 and 101 
("'loiQO^ ), 127 (I - ), and 191 (I0 4 - ) were monitored. 

The mass spectra recorded reflect the composition of ions 
that belong to the interface layers of the ESI jet at the time it is 
aerodynamically sheared by the nebulizer gas. Previous related 
studies characterized the outermost layers of aerosolized 
droplets using a grounded capillary (electrode) ESI 
probe. 18,24 However, it was shown that the behavior of 
conducting and nonconducting electrospray capillaries pro- 
duced similar results. 32 The experimental results presented in 
the next section are validated by showing that the abundance of 
seven anions, as reflected by their MS ion count, quantified 
from equimolar solutions follows a normal Hofmeister 
sequence. 18,24 Therefore, the enrichment is not similar for 
each ion, as expected at the air— water interface, and the 
fractionation pattern is specifically affected by the addition of 
surfactants. 

The flow of a solvent pump (Lab Alliance) at 150 /tL/min, 
stabilized with a back pressure regulator, was mixed through a T 
connector with a 50 ^L/min flow of the sample in a syringe 
pump (Harvard Apparatus, Elite Series) and directed to the 
ESI-MS (Thermo Scientific, MSQ_ Plus). The solvent back- 
ground was subtracted from the mass spectrum of the sample 
acquired at fixed time intervals of at least two minutes. 
Reported data are the average of duplicate experiments. 

Present experiments are also supported by confirming the 
lack of any kind of correlation between the current applied to 


the ESI capillary and the fractionation process observed. 
According to Kohlraush’s law, each ion in the 100 nM 
equimolar (dilute) solution carries its portion of the total 
conductivity without being affected by any of the other ions. 40 
No dependence is observed between the ionic limiting 
equivalent conductivity for each ion (A°) and fractionation 
(see Supporting Information, Figure S1A). Therefore, under 
the conditions chosen to study the fractionation process, the 
major controlling factors are the fast mixing between the 
nebulizing gas and the liquid (t'gas/Hiquid = 119 X 10 3 ) during 
pneumatically assisted aerosolization and the Coulomb fission 
of droplets showed in the results section. In addition, the same 
fractionation behavior for all the ions was observed for the 
nebulizer voltage range between 0.5 and 2.0 kV (Supporting 
Information, Figure SIB). 

A solution composed of 1.3 mM Nal, 0.05 mM KI, and 0.02 
mM Csl, in 50:50 water— 2-propanol (100%, Fisher Optima) 
solvent, was used to calibrate the mass spectrometer within the 
mass range m/z 22.9898—1971.6149 amu. The calibration 
conditions were drying gas temperature, 350 °C; nebulizer 
voltage, 3 kV; cone voltage, 75 V; and nebulizer pressure, 70 
psi. 

■ RESULTS AND DISCUSSION 

Experiments with Equimolar Solutions. The mass 
spectrum of a 100 nM equimolar solution of NaCl, NaN0 2 , 
NaBr, NaN0 3 , Nal, NaC10 4 , and NaI0 4 in water is shown in 
Figure 1A, relative to the ion intensity of I0 4 - . The negative 
ions present in the solution do not display uniformity in the ion 
counts detected due to the process of interfacial fractionation. 
Fractionation enhances the concentration of I0 4 - (the most 
abundant ion), C10 4 - , and I - over N0 3 - and Br“. 
Simultaneously, fractionation practically suppresses the signals 
of N0 2 _ and Cl - , only observed in the inset of Figure 1. The 
mass spectrum of Figure 1 can be normalized as an anion 
enrichment or interfacial fractionation, f x , calculated from 

j- ^m/ z,X~ 

2 Im/z,X~ ( 1 ) 

where I m/z is the sum of the ion count for all the isotopes of 
each anion X - , e.g., (I 79 + I 81 ) for Br - , and the denominator is 
the total ion count (Table l). Although the individual ions 
show a linear dependence of I m/z on concentration, f x is 
independent of concentration between 50 and 200 nM, because 
it is a relative measurement that reflects the composition of the 
outermost layers of the droplets instead of the composition of 
the bulk solution. 18,24 

In order to classify the dominant mechanism for the selective 
enrichment and depression of anions in the interface, it is useful 
to first explore how the different ion properties affect the 
process. Anion size has been proposed to control f x ~ in related 
studies 18, “ 4 that explored its dependence on ion radius ( R x ~ ) to 
the first power. The implication is that f x ~ should depend on 
the free energy change associated with the rejection of anions 
from the liquid bulk to the air— water interface through a 
(negative) exponential function, f x ~ cc exp[— Bulk ^ Interface AG/ 
kT], 18 that predicts a linear correlation between ln/ x - and R x -. 
However, the possibility of several anions interacting at the 
interface has been dismissed so far. From the viewpoint of 
anions repelling each other, the interaction energy should be 
inversely proportional to the square of the ion radius, and a 
straight line with the best overall correlation can be expected in 
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Figure 1. (A) Spectra of ESI-MS of 100 nM solution of NaCl, NaBr, 
Nal, NaNOj, NaN0 3 , NaC10 4 , and NaI0 4 at pH 6. Ion count values 
normalized relative to I0 4 - count are reported in Table 1. (B) Anion 
fractionation factor, f x -, versus reciprocal square of anion radius (from 
ref 42), R y -I . The correlation coefficient (r 2 ) of the linear fit shown 
with dashed lines are (black open circle) 0.948 for H 2 0, (red open 
triangle) 0.948 for D 2 0, (green open diamond) 0.923 for CH 3 OH, 
and (blue open square) 0.942 for CH 3 CH 2 OH. 


a plot of ln/ x - vs R x -~ 2 . Herein, we show that our experiments 
are better correlated with R x ~~ 2 , what indicates that the 
enrichment of anions at the interface also depends on the 
Coulomb interactions among ions. When considering ion 
polarizability (a), a direct increase of f x is expected. 41 In the 
case of anion dehydration free energy (AG dehyd ), 18,28 (Table l), 
and given the inverse proportionality AG dehyd ex R x - -1 , 42 the 
best correlation expected should linearly link a decrease of In 
fx with AG dehyd . 

The fact that the largest and more polarizable I0 4 - (R I0 ~ = 
249 pm; a = 7.69 X 10 -30 m 3 ) 42 is the most abundant ion at the 
interface suggests a correlation between enrichment with anion 
radius (R x ~) and polarizability (a). This observation agrees with 


the smallest and less polarizable Cl - (Ref = 181 pm; a = 3.42 
X 10 -30 m 3 ) 42 being the less enriched ion at the interface or, in 
other words, the most enriched ion in the solvent bulk. 
However, the Coulombic repulsion among anions is the 
controlling factor during bubble fission. Figure IB shows the 
linear relationship of In f x ~ vs the square of the reciprocal of 
anion radius ( R x ~~ 2 ) for solutions in pure H 2 0, D 2 0, CH 3 OH, 
and CH 3 CH 2 OH, respectively. The correlation of In f x ~ with 
R x - (0.90) is smaller than with R x ~~ 2 (0.95) or R x ~~ l (0.93) 
(Figures 1 and S2, Supporting Information). Indeed, the best 
fitting corresponds to R x ~~ 2 . These results reflect that the 
fractionation process observed is due to the Coulomb fission of 
microdroplets with excess negative charge and many-anion 
repulsion. The correlation of In f x with R x ~ and R x ~ , 
indicate a dependence on Coulomb’s energy and dehydration 
free energy (to be discussed below), respectively. It is apparent 
that the ratios of the linear regression slopes in Figure IB, 
m H2o/ m D 2 o = —4.46 X 10 5 pm -2 /— 4.49 X 10 s pm -2 = 0.993 
and WcH 3 OH/ m CH 3 CH 2 OH = — 4.26 X 10 s pm "/— 4.27 X 10 3 
pm -2 = 0.998, as well as the correlation coefficients ratios 
r H 2 o 2 / r D 2 o 2 = 0.948/0.948 = 1.00 and t"cH 2 OH 2 / r CH 3 CH 2 OH 2 = 
0.923/0.942 = 0.980, are practically identical within exper- 
imental error. The previous observations also indicate the lack 
of any isotopic effect when substituting H 2 0 for D 2 0, which is 
a demonstration that we are dealing with a pure physical 
fractionation process, where no chemical reaction is taking 
place, as expected in solutions of these electrolytes. 

Polarizability. Figure 2 shows the In f y vs polarizability of 
the ions (a x ~), which reflects the relationship between 
fractionation and the ability of the electronic shells of the 
ions to undergo deformation in an electric field. 42 The poor 
correlation coefficient between ln/ x - and a x (r H2 o 2 = 0.679) is 
closer to the correlation between a x and R x (r aXx _ 2 = 0.592) 
also shown in Figure 2 but considerably smaller than those 
observed with R x ~ 2 in Figure IB. Indeed, the minor 
contribution of polarizability to the fractionation process can 
be explained by its own dependence on ion radius. 

Effect of Surfactant. The effect of polysorbate 20 on anion 
fractionation, as nonionic surfactant, seems negligible for 
surfactant concentrations below 10 -5 mol L -1 (Figure 3). 
The addition of the surfactant below the critical micelle 
concentration (CMC) 43 to the electrolyte solution decreases 
the total ion count 12 (the ion count for Cl - vanishes); the air- 
water interface is altered, modifying its interfacial fractionation 
at even low bulk concentrations. The surfactant occupies the 
interface in an oriented fashion and forms micelles in the bulk 
phase at the CMC 43 = 8.0 X 10 -s M of polysorbate 20. 44 Figure 


Table 1. Ion Properties and Fractionation of Anions at the Air— Water Interface" 


X - 

R x - (pm) 

a (10 -30 m 3 ) 

AG dehyd (kj moF 1 ) 

2° (10~ 4 m 2 S mol -1 ) 

I x - C (equimolar) 

f x ~ c (equimolar) 

/x“//cf d (seawater) 

cr 

181 

3.42 

347 

76.4 

1.19 X 10 -3 

6.36 X 10 -4 

1.00 

no 2 ~ 

192 

3.45 

339 

71.8 

6.65 X 10 -3 

3.54 X 10 -3 

1.93 X 10 -4 

Br - 

196 

4.85 

321 

78.1 

4.39 X 10 -2 

2.34 X 10 -2 

9.12 X 10 -2 

no 3 - 

20 6 b 

4.13 

306 

71.5 

3.63 X 10 -2 

1.93 X 10 -2 

9.40 X 10 -4 

r 

220 

7.51 

283 

76.4 

0.240 

0.128 

2.96 X 10 -4 

cio 4 - 

240 

5.06 

214 

67.4 

0.549 

0.292 

/ 

io 4 - 

249 

7.69 

n.a. e 

54.6 

1.00 

0.533 

1.37 X 10 -3 


“Experiments were carried out in H 2 0. Key: X - , anion; R y , radius; a, polarizability; A G dt .| ly j, dehydration free energy; 2°, limiting equivalent ionic 
conductivity; I Xl normalized ion count; f x , fractionation factor. R x -, a, AG dehyd , and 2° values are from ref 42. “This value corresponds to the 
equatorial radius. ‘'Equimolar concentration data for H 2 0 in Figure 1. ^Seawater molar ratio, ‘n.a., data not available, fflelow the detection limit. 
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Figure 2. Anion polarizability (from ref 42) versus (black open circle) anion fractionation factors, / x , and (blue open triangle) anion radius, R x . 
The linear fit yields correlation coefficients (r 2 ) of 0.679 for f x ~ and 0.S92 for R x . 

(depressed) in the total number of ions still present at the 
interface. The opposite trend of C10 4 - and I0 4 - is explained 
below. 

The effect of surfactant on f x ~ can be explained in terms of 
salting out or salting in relative to the hydrophobic groups of 
polysorbate 20 in water. 43 Since the hydrophilic groups of the 
surfactant molecules are in contact with the aqueous phase in 
both the monomeric and micellar forms of the surfactant, while 
the hydrophobic groups are in contact with the aqueous phase 
only in the monomeric form, the fractionation of anions is 
affected by the addition of surfactant. The effects that 
electrolytes experience is close to be canceled out in the 
monomeric form below [polysorbate 20] = 10 -5 M. Anions 
with low polarizability ( a < 5.1 X 10 -3 ° m 3 , Table l), such as 
N0 2 - , N0 3 - , Br - , and C10 4 - are water structure makers. They 
salt out the hydrophobic groups of the surfactant close to the 
CMC and are enriched compared to those that salt in. Anions 
with large polarizability (a > 7.5 X 10 -3 ° m 3 ), I - and I0 4 - , are 
water structure breakers; they salt in the hydrophobic groups of 
the surfactant and are depleted above the CMC. Thus, the 
enrichment ratio of anions with and without surfactant, / X "(S)/ 
fx~( 0), above the CMC is affected by the polarizability and 
geometry of anions. 43 

Effect of Alcohols. Figure 4 shows that f x ~ measured in 
water— methanol and water— ethanol mixtures is only weakly 
dependent on the composition of the solvent for molar 
fractions 0 < x HlQ < 1 (x r _oh = 1 — 2C H2 o)- Sum-frequency 
generation vibrational spectroscopy measurements have shown 
that, for increasing 2 C CH3 oh — 0-2 , the dangling —OH of water 
gradually disappears. 35 A similar behavior is expected for 
ethanol. 36 The configuration includes the polar —OH group of 
the alcohol pointing into the bulk phase, satisfying hydrogen 
bonding and minimizing dipole repulsion. Independently of the 
surface structure tested in our experiments by gradually 
substituting the interfacial hydroxyl groups of water by methyl 
and ethyl groups, 34-37 anions selectively occupy the interface. 
Indeed, all anions should be oriented to minimize the distortion 
of the interface. Anions can repel each other and distribute in a 
way that provides the most gradual possible transition from the 
bulk solution to air. In consequence, the molecules of the 



Figure 3. Ratios of anion fractionation factors in the presence and 
absence of surfactants, / x (S)// x (0), at variable concentration of 
polysorbate 20: (blue x) N0 2 - , (green filled square) NO, - , (red filled 
upside-down triangle) Br - , (black filled diamond) C10 4 - , (purple 
open triangle) I - , and (dark green filled circle) I0 4 - . 


3 shows that, at low surfactant concentration, / C io 4 “) / io 4 ~; an d 
/ N o 2 - are close to their values in water, f(0) ao ~, and 

/(0)no 2 T while /p and/ Br ~ are depressed. Considerable changes 
are only observed for [polysorbate 20] > 10 -s M. At higher 
polysorbate 20 concentration, there is an interfacial enrichment 
in the population of N0 3 - > N0 2 - > Br - > C10 4 - in the 
surface layer and a fast decrease of the largest and more 
polarizable I0 4 - that together with I - disturb the stability of the 
surface active layer. In the presence of 1 mM surfactant, the 
strong correlation of In f x with R x ~ 2 drops to r 2 = 0.580. 
Although the total population of ions (or its density) occupying 
the air— water interface region decreases in the presence of 
surfactants, the fractionation of anions that have a planar 
molecular geometry f NOj ~ and /no 2 “ is enhanced at the air- 
water interface with surfactant as compared to pure water 
/no 3 “(0) and/ NO2 -(0). Other ions that show opposite behaviors 
at [polysorbate 20] = 10 -J M are C10 4 - (enhanced) and I0 4 - 
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Figure 4. Anion fractionation factors, f x , for 100 nM solution of 
sodium salts (top panel), interface thickness (from ref 37), and 
dielectric permittivity of the solvent (from ref 45), £, (bottom panel) 
at variable (A) methanol and (B) ethanol molar fractions: (blue x) 
N0 2 “, (green filled square) N0 3 - , (red filled upside-down triangle) 
Br - , (black filled diamond) C10 4 - , (purple open triangle) I - , (dark 
green filled circle) I0 4 - , and (blue-green filled circle) Cl - . 


solvent at the interface should still be oriented in a way so that 
their mutual interaction energy will be a maximum. 34 

For the methanol hydrophobic interface (as x CH 0H — » 1.0), 
when the surface number density reaches ~3 methanol 
molecules nm -2 , 35 there is no possibility of water molecules 
interacting with anions, and the governing intermolecular forces 
involve methanol hydrogen bonding and ion-dipole inter- 
actions. The change in the interfacial layer thickness or thermal 
roughness with bulk concentration, evaluated as a function of 
the interface surface tension reported by Sung et al., 37 is plotted 
in the bottom panels of Figure 4. Both water— methanol and 
water— ethanol interfaces have a strongly oriented outer layer, 
with a hydrophobic — CH 3 and — CH 2 — CH 3 packing structure 
of thinner width than the radius of the methanol and ethanol 
molecules, respectively. 34 ’’ 6 The fractionation factor f x ~ barely 
changes when the thickness of the outermost layer of the 
interface 37 is greater than ~5 A with increasing alcohol bulk 
mole fraction larger than 0.2. 

At the interface, a solid spherical ion of radius R x ~ and charge 
q may have one of its hemispheres submerged in the liquid bulk 
(a continuous dielectric medium) to a depth z and the other 
exposed to air. 20 The effective dielectric permittivity of the 
interfacial layer should gradually decrease from that of the 
hydrogen bonded solvent to air. Figure 4 (bottom panels) also 
presents the decrease in the dielectric permittivity of methanol 
and ethanol— water mixtures for increasing bulk alcohol molar 
fraction. 45 When considering the interfacial liquid slab, 
electrolytes strongly interact with the surrounding solvent 
molecules through ion— dipole interactions and occupy a 
position to decrease the dielectric permittivity of the solvent 


^solvent- 46 Therefore, the dielectric conductivity of the hydrated 
ions, e x ~, lays between the values for the interfacial layer and 
the solvent e(z)soivcm C < e x - < Solvent, where t' ) solvent * t?ai r 
1 asz-> 0. 46 A competing factor to keep the ions in the bulk is 
their large dehydration free energy A G dehyd) 47 although its 
contribution seems minor (q n / x _, AG Jlh , d 2 = 0-786) based on the 
data shown in Figure S3A, Supporting Information. Atomistic 
simulations of ion solvation in water, dependent on the 
enthalpic component, suggest the surface preference of halides 
based on their energy minima at the surface. A partially 
desolvated anion favors the water— water interactions more 
than in the fully solvated state. 28 In addition, anion 
fractionation appears to be independent in the range of pH 
6.8— 8.6, a relevant interval for seawater, and does not change 
when substituting the countercation Na + for Li + , K + , or Cs + in 
100 nM equimolar solutions containing Cl - , Br - , and I - 
(Figure S4, Supporting Information). 

Seawater Mimic Solutions of Monovalent Anions. In 
spite of the observed enrichment of iodide for the air— water 
interface of equimolar solutions, given the large concentration 
difference with actual seawater composition, two mimic 
samples of the monoanions were studied. For the first sample 
I, the concentration of sodium salts were 100 fiM Cl - , 0.44 nM 
N0 2 - , 1.8 nM N0 3 - , 121 nM Br - , 41 pM I - , 0.31 pM C10 4 - , 
and 100 nM I0 4 - . All concentrations are scaled down to a 
reference of 100 /.iM Cl - , and the molar ratios of monovalent 
anions in surface seawater Cl - /N0 2 - /N0 3 - /Br - /I - /C10 4 - / 
I0 4 - were kept fixed to 100:4.4 X 10 -4 :1.8 X 10 -3 :0. 12:4.1 X 
10 -s :3.1 X 10 -7 :0.10. 29-33 The aerosolization of this sample 
yields f x values (normalized to f a ~ = l) for seawater molar 
ratios of anions (Figure S5A, Supporting Information). The 
experimental result shows that for actual seawater, the ratio y = 
fa~/f Bf ~ 26 indicates that Cl - is 26 times more enriched than 
Br - due to its larger concentration in solution. The second 
seawater mimic sample II keeps the molar ratios for a 10-times 
dilution of anions in seawater, with the exception of Cl - that is 
1000 times diluted: 545 Cl - , 239 nM N0 2 - , 957 nM 
N0 3 - , 66.24 /(M Br - , 22.2 nM I - , 0.17 nM C10 4 - , and 100 nM 
I0 4 - . 29 33 This sample was useful to confirm that Br - becomes 
the second most enriched species (Figure S5B, Supporting 
Information). The molar ratios of anions C1 - /N0 2 - /N0 3 - / 
Br - /I - /C10 4 - are fixed to 100:4.4 X 10 -2 :0.18: 12.15:4.1 X 
10 -3 :3.1 X 10 -5 . Fractionation factors in samples I and II vary 
according to the change in relative concentrations. 

Figure 5 shows a more realistic scenario of seawater 
fractionation obtained by multiplying f a ~ in sample II by the 
scaling factor y obtained from sample I, where f C f = 1 means 
that Cl - is by far the most enriched ion transferred to the 
atmosphere during bubble bursting, and/ Br - = 0.091 accounts 
for Br - being the second most enriched (9.1% of Cl - ) anion 
(Table l). The bulk molar ratio [Br — ] / [l — ] = 132 is 2.3 times 
smaller than the relative excess of Br - at the interface fsf/ff = 
9.10 X 10 -2 /2.9 6 X 10 -4 = 308, confirming that the interfacial 
enrichment of I - for equimolar solutions is not followed at 
relevant seawater molar ratios. Indeed, the role of concentration 
has to be considered in any attempt to explain the enrichment 
of ions at the air— water interface. 

■ CONCLUSIONS 

This study addresses the fundamental problem of the selective 
transfer of anions occurring during aerosol spray formation. 
The enrichment of anions at the air— water interface 
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Figure 5. Anion fractionation factor, f x ~, normalized relative to / cl -, for 
molar ratios of anions found in seawater. 
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represented by the fractionation factor shows an increasing 
linear correlation with ion size that follows the order R x ~ < 
R x -_1 < Rx~ 2 . An explanation to the latter pattern is that the 
surface tension of the droplets is overcome by a series of many- 
body interactions among the solvated anions. These many-body 
interactions are stronger than the work of adhesion among 
solvent— solvent molecules with their associated pairwise 
interaction energies. The series of cascade Coulombic fission 
that droplets undergo, and occur a posteriori of the surface 
sampling stage, produce a pure physical enrichment of larger 
and less repulsive anions at the interface. Although the ionic 
limiting equivalent conductivity is not correlated with the 
selective enrichment of anions at the interface, minor 
collinearity contributions from other ion properties such as 
polarizability and dehydration free energy exist. The minor 
correlation of anion enrichment factors with polarizability can 
generally be explained based on the dependence of polar- 
izability with anion size (r aRx _ 2 = 0.592). In addition, 
polarizability is not a significant predictor to be included in 
the model. The correlations between R x ~ and AG dehyd 
( r AG de h d,«x— 2 = 0.960) indicates that the combination of these 
ion properties to explain the trend of In f x ~ represents 
redundant information and little can be gained by including 
AG dehyd in the model. For equimolar solutions, anion 
enrichment factors at the air— water, air— methanol, and air- 
ethanol interfaces are mainly dependent on R x ~~ 2 - The 
depletion of anion species in the presence of a nonionic 
surfactant above the CMC, as [X“] inlerface — » 0, follows a 
polarizability trend. The experimental result for seawater 
mimics of monovalent anions distinguishes the effect of 
concentration on anion enrichment factors at the air— water 
interface. For the latter sample, the enrichment factor ratio 
if c f // Br = 1 1) is 76 times smaller than the concentration ratio 
of these species in seawater ([Cl _ ]/[Br“] = 833). Indeed, the 
fractionation trend observed is somewhat biased by its 
dependence on size as R x ~~ 2 . Contrarily to the results for 
equimolar solutions, when the seawater molar ratios are fixed in 
the bulk solution, the aerosolized samples show that the air- 
water interface is 2.3 times richer in bromide than iodide. 
Present results reveal the relative position that anions occupy at 
the air— water interface and provide direct insight into the 
source of halides for the production of tropospheric RHS that 


participate in bromide explosion events in the early Arctic 
spring. 
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